Astrophysics and Space Science 

DOI 10.1007/sXXXXX-XXX-XXXX-X 



The.Effe.ct of R.adiation Pressure on the 
Equilibrium Points in the Generahzed 

Photogravitational Restricted Three Body 

Problem 

Badam Singh Kushvah-*^ 



© Springer- Verlag ' 



Abstract The existence of equilibrium points and the 
effect of radiation pressure have been discussed numeri- 
cally. The problem is generalized by considering bigger 
primary as a source of radiation and small primary as an 
oblate spheroid. We have also discussed the Poynting- 
Robertson(P-R) effect which is caused due to radiation 
pressure. It is found that the collinear points Li, L2, 
deviate from the axis joining the two primaries, while 
the triangular points L4, L5 are not symmetrical due 
to radiation pressure. We have seen that Li, L2, are 
linearly unstable while L4, are conditionally stable 
in the sense of Lyapunov when P-R effect is not consid- 
ered. We have found that the effect of radiation pres- 
sure reduces the linear stability zones while P-R effect 
induces an instability in the sense of Lyapunov. 

Keywords Radiation Pressure: Equilibrium Points: 
Generalized Photogravitational:RTBP:Linear Stability: 
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1 Introduction 

Three-Body problem is a continuous source of study, 
since the disco very of its non-integrability due to 
Poincard (jl892l ). Many of the best minds in Math- 
ematics and Physics worked on this problem in the 
last century. Regular and chaotic motions have been 
widely investigated with any kind of tools, from ana- 
lytical results to numerical explorations. The restricted 
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three body model(RTBP) is insoluble, although spe- 
cific solutions exist, like the ones in which the space- 
craft is positioned at one of the Lagrangian points. 
The five singular points of Jacobian function are called 
equilibrium points, Lagrangian points and most fre- 
quently are referred to as the equidistant (triangular) 
and collinear (straight line) solutions. There are many 
periodic orbits in the restricted three body problem. 
One of the most famous, discovered by Lagrange, is 
formed by an equilateral triangle. In "Earth-Moon- 
Space Station" model, if the Moon is not too massive, 
this orbit is thought to be stable. If the space station 
is pushed a bit to one side (in position or velocity ), 
it's supposed to make small oscillations around this or- 
bit. The two kinds of triangular points are called L4 
and L5 points. We study the motion of three finite 
bodies in the three body problem. The problem is re- 
stricted in the sense that one of the masses is taken to 
be so small that the gravitational effect on the other 
masses by third mass is negligible. The smaller body 
is known as infinitesimal mass and remaining two finite 
massive bodies as primaries. The classical restricted 
three body problem is generalized to include the force 
of radiation pressure, oblateness effect and Poynting- 
Robertson(P-R) effect. The solar radiation pressure 
force Fp is exactly opposite to the gravitational attrac- 
tion force Fg and change with the distance by the same 
law it is possible to consider that the result of action 
of this force will lead to reducing the effective mass of 
the Sun or particle. It is acceptable to speak about a 
reduced mass of the particle as the effect of reducing 
its mass depe n ds on the p roperties of t he par ticle itself. 
Chernikovl (jlQTdh and ISchuermai] (|l980[ ) discussed 



the position as well as the stability of the Lagrangian 
equilibrium points wh en radiat i on pre ssure, P-R drag 



force are included. iMurravl (|l994r ) systematically 



discussed the dynamical effect of general drag in 
the planar circular restricted three body problem. 



2 



Ishwar and Kus livah ( 2006 ) examined the linear stabil- 
ity of triangular equilibrium points in the generalized 
photogravitational restricted three body problem with 
Poynting-Robertson drag, and L5 points became 
unstable due to P-R drag which is very remarkable 
and important, where as they are linearly stable in 
classical problem when < m < URny .th = 0.03852. 
Kushvah. Sharma. and Ishwan ( 2007al lbl)d) examined 



normalization of Hamiltonian they have also studied 
the nonlinear stability of triangular equilibrium points 
in the generalized photogravitational restricted three 
body problem with Poynting-Robertson drag, they have 
found that the triangular points are stable in the non- 
linear sense except three critical mass ratios at which 
KAM theorem fails. 



2 Equations of Motion 



Povntinj ( 1903 ) has stated that the particle such as 
small meteors or cosmic dust are comparably affected 
by gravitational and light radiation force, as they ap- 
proach luminous celestial bodies. He also suggested 
that infinitesimal body in solar orbit suffers a gradual 
loss of angular momentum and ultimately spiral into 
the Sun. In a system of coordinates where the Sun is 
at rest, radiation scattered by infinitesimal mass in the 
direction of motion suffers a blue shift and in the op- 
posite direction it is red shifted. This gives rise to net 
drag force which opposes the direction of motion. The 
proper rela tivistic treatment of this problem was for- 
mulated by I Robert soni (|l937l ) who showed that to first 
order in — the radiation pressure force is given by 



V.RR V 

c 



(1) 



Where Fj,= — denotes the measure of the radi- 
ation pressure force, R the position vector of P with 
respect to radiation source Sun 5, V the correspond- 
ing velocity vector and c the velocity of light. In the 
expression of i^p, L is luminosity of the radiating body, 
while TO, p and s are the mass, density and cross section 
of the particle respectively. 

The first term in equation ([T]) expresses the radiation 
pressure. The second term represents the Doppler shift 
of the incident radiation and the third term is due to the 
absorption and subsequent re-emission of the incident 
radiation. These last two terms taken together are the 
Poynting-Robertson effect. The Poynting-Robertson 
effect will operate to sweep small particles of the solar 
system into the Sun at cosmically rapid rate. We con- 
sider the barycentric rotating co-ordinate system Oxyz 



relative to inertial system with angular velocity uj and 
common z-axis. We have taken line joining the pri- 
maries as a;-axis. Let toi,TO2 be the masses of bigger 
primary(Sun) and smaller primary(Earth) respectively. 
Let Ox, Oy in the equatorial plane of smaller primary 
and Oz coinciding with the polar axis of 7712. Let re, Vp 
be the equatorial and polar radii of TO2 respectively, r be 
the distance between primaries. Let infinitesimal mass 
TO be placed at the point P{x, y, 0). We take units such 
that sum of the masses and distance between primaries 
is unity, the unit of time i.e. time period of toi about 
TO2 consists of 27r units such that the Gaussian con- 
stant of gravitational = 1. Then perturbed mean 
motion n of the primaries is given by = 1 -|- 



3A2 



where A2 



5^2 



is oblateness coefficient of TO2. Let 



then 1 — /i 



with TOi > TO2, where 



p, is mass parameter. Then coordinates of toi and TO2 
are (— /i, 0) and (1 — /i, 0) respectively. Further, in our 
consideration, the velocity of light needs to be dimen- 
sionless, too, so consider the dimcnsionless velocity of 
light as Cd = c which depends on the physical masses of 
the two primaries and the distance between them. In 
this paper, we set = 299792458,/^ = 0.00003 for 
all numerical results. 

In the above mentioned reference system the total 
acceleration on the particle P is as follows 

a' = a -\- 2lu X V -\- uj X [uj x f) 
_ (l-/^)rl nf"2 3 ^^2^1 



+ - 



(1 - ^)(1 - gi) frl {ri + LO X fl).flfl 



Cdr'i 



ri -\- Q X fi 

Cd 



(2) 



where 



fi = xi + yj , V = xi + ijj, a = xi + yj , uj = nk, 
fi = {x + p)i + yj, fi = {x + p, - l)i + yj, 
rl = {x + fif + 2/2, ri = (x + M - 1)' + y^ 

UJ X V = —n[yi — xj), uj x (uj x f) = 71^ [xi -\- yj), 
r'l +nk X fi = {x- ny)i + [y + n{x + ^)] j, 
(^fl +nkx fi^ .fi = [x[x + ^) + yy] 

Substituting all these values in above relation ( [21) we 
get as follows: 
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(i — 2ny — n^x) ? + (y + 2nx — rP"y) j 

_ (1 - /i)fi _ ^JT2 _ 3^f2^2 ^ (1 - M)gl^l 

(1 - - qi) / ( ri + X fl).flfl 



- uj X ri 



Cd 



(l-/i)giri ^lr2 3^f2A2 



is the Jacobi's constant. The zero velocity curves[see 
([2l),([3])] are given by: 



Ci = 2U{xi,yi) 



(6) 



SufRx i = 1,2,3,4,5 correspond to respective i*'' La- 
grangian equihbrium point Li. Using above relation we 
have determined Ci « 3.02978, Ca ~ 4.04133, C3 « 
3.53607. The values of C4(« C5) are shown by dif- 
ferent curves (1 — 6) in figure { [IJ and table ( [I]) for 
various values of qi , A2 



Wi ( (r'l + a; X ■f{).fifi r 1 + x rl 



Cdrf 



1 



^1 ~ — — ' substituting the values of fi , r2 and 
comparing the components of i and j we get the equa- 
tions of motion of the infinitesimal mass particle in xy- 
plane. 



X — 2ny 
y -\- 2nx 



(3) 

(4) 



where 



2 _ (l-/i)gi(a^ + A^) _ /^x-H/^- 1) 



3 /iA2(a: -|- ^ - 1) 
"2 Vl 

^1 I ^1 
_ 2 (l-/^)giy My 3 ^A2y 
rj* r2 2 

~—T\\[i^ + ^')^ + yy\ + y + n{x + 
^1 L ''i J 



where 



{x^+y'^) {l~^i)qi u llA 



+ 



+ — + 



r2 



[x + + yy 
+Wi < ^ n arctan 



x + ^ 



5) 



qi = 1 — -gf- is a mass reduction factor expressed 
in terms of the particle radius a, density p radia- 
tion pressure efficiency factor x (in C.G.S. system): 
gi = 1 — ^'^^^p — X- The assumption qi = constant 
is equivalent to neglecting fluctuations in the beam of 
solar radiation and the effect of the planets shadow, 
obviously f/i < 1 . The energy integral of the problem 
is given by C = 2U — x^ — y^, where the quantity C 



3 Existence of Equilibrium Points 

What are they? What are "Lagrange points", also 
known as "libration points" or "L-points" or "Equilib- 
rium Points"? These are all jargon for places where 
a light third body can sit "motionless" relative to two 
heavier bodies that are orbiting each other, thanks to 
the force of gravity. The unstable Lagrange points - 
labeled Li, L2 and L3 -lie along the line connecting the 
two large masses. The conditionally linearly stable La- 
grange points in classical case are labeled L4 and L5 
-form the apex of two equilateral triangles as in figure 
([3]), that have the large masses at their vertices. 

3.1 Collinear Equilibrium Points 

To investigate the Equilibrium Points divid the orbital 
plane Oxy into three parts with respect to the primaries 
X < — /i, 1 — ^ < X and — /i < x < 1 — fi, Ux = Uy = 




Fig. 1. 

(2)^2 
(6)^2 



4.5 



- The Jacobi's constant C4 — qi for (1)^2 = 

0. 2 {3)A2 = 0.4 (4)^2 = 0.6 (5)^2 = 0.8 

1, 0<9i < 1 &Ai = 0.00003 
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Table 1 Jacobi's constant C at ^4 


A2 


C4 : 91 = 1 


C4 : gi = 0.75 


C4 : 91 = 0.5 


C4 : 91 = 0.25 


C4 : 91 = 


0.0 


2.99997 


2.47643 


1.88988 


1.19058 


Indeterminate 


0.2 


3.27697 


2.70632 


2.06659 


1.30303 


Indeterminate 


0.4 


3.51014 


2.89885 


2.21377 


1.3961 


Indeterminate 


0.6 


3.71975 


3.0698 


2.34214 


1.47515 


Indeterminate 


0.8 


3.98328 


3.30328 


2.53573 


1.61041 


Indeterminate 


1 


4.62593 


4.06149 


3.42503 


2.62291 


Indeterminate 




then from equations ( [3]) and ( 
^1-^)91 ^l 



2 



= 0. 

' nWi 



rl + -A2 = 0, 



rl [(1 - y)yn'^ + nWi] - n{x + n)Wiri = (1 - ^J)qly 



Fig. 2.— CYsx-y when A2 = 0.0024, < 91 < 1 & 
^ = 0.00003 



0.5 



-0.5 




Fig. 3. — Contour plot shows the Jacobi's constant C 
Vs x-y when A2 = 0.0024, < 91 < 1 & /x = 0.00003 



r-2 



1/3 



1 - 



nWi 



My 2 



-1/3 



(7) 
(8) 



From above, we obtained: 



'2L 

r> 2 



M± 

2/3 



1 + 



nWi 



I- ^l± 



2{l~y)y 



+ 



3^2 



-1 -2/3 



1 



-2/3 



1/2 

(9) 

1 1/2 

(10) 



Using equations ( [H [TO]) the position of Li,L2,L3 
is presented graphically in figures ( [H [5]) when A2 = 
0.0024, 0<9i<l&M = 0.00003. We have seen that 
the equilibrium points are no-longer collinear with the 
primaries. The values of y are positive. We observe 
that the coordinates of L3 are the functions of 91 and 
A2, corresponding different curves labeled by (1 — 4) 
are presented as in figures ( [5|). Now when 1 — /i < x, 
there exists an equilibrium point L2 this can be seen in 
figures (ElE]). This equilibrium point is also found away 
from the Ox axis. It is clear from the figure that L2 is 
a function of 91 and ^2- When — ^ < a; < 1 — /i, there 
are almost three equilibrium points they are Li, L4 and 



The Effect of Radiation Pressure on the Equihbrium Points in the GeneraUzed Photogravitational Restricted Three Body Problem 5 













/ 0.5 














i Li 


lE 








-0.5 


0.5 




1.5 















0.01 
0.0075 
0.005 
0.0025 


-0.0025 
-0.005 
-0.0075 



■ 1 1 1 1 1 1 1 1 1 


u 


Cj =1*413332 


u 

1-^ 
".■1 


L2 


■i ■ \ : 


".A 

"."a. 
Vj 




•wi 
H 





■1.9 



1.95 



2 

X 



2.05 



2.1 



Fig. 4. — Figure views the Lagrangian equilibrium 
points in classical case(when A2 = 0,qi = = 
0.00003) 



Fig. 7. — Figure views the magnified region of L2, when 
A2 = 0.0024 & Ai = 0.00003, < gi < 1 
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Fig. 5. — The position of Li, L2, L3, L4 and for 
= 1 (2)gi = 0.75 (3)gi = 0.25 {A)qi = 0, 
(5)0 < qi < 1 when A2 = 0.0024 & ^ = 0.00003 
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Fig. 8. — Position of the Lagrangian equilibrium points 
L3 when A2 = 0.0024 & /i = 0.00003 for (l)gi = 1 
{2)qi = 0.75 (3)gi = 0.25 (4)gi = 




X 



Fig. 6. — Position of the Lagrangian equilibrium point 
L2 when A2 = 0.0024 k fi = 0.00003, < 91 < 1 



Fig. 9. — Position of the Lagrangian equilibrium points 
Li,L4,L5 when A2 = 0.0024 & ^ = 0.00003, for 
(l)gi = 1 (2)gi = 0.75 (3)gi = 0.25 (4)gi = 
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Ls as in figure ([9]). It is clear from the figure, that the 
Li is not colhnear any more. L4 is positioned above 
the Ox axis while L5 h e s below it. All these results are 
simila r withlSzebehelvl (119671). iRagos and Zafiropoulos 
( 19951 ). Papadakis and Kanavos ( 2007 ) and others. 



3.2 Triangular Equilibrimn Points 

For the triangular equilibrium points y ^ 0, Ux = Uy = 
0, then from equations ([3|) and (|4|) we get as follows: 



ri^x 



3fiA2{x + ^-l) Winy 



+ 



0, 



n'^y 



Wi 



'2 '1 

fiy 3fiA2y 



.,{x + ^) = 0. 



(11) 



(12) 



0.4 



0.2 



[x + ^) + 



nWi 




^< n 



1 + 



3A. 



^1 - 



[1 + 



3A2 



V 



nWi 



0. 



From equations ( fTTj) and ( [T2|) . we obtained as follows: 

nWi. (13) 



-^(1 + 1%' 



My 



In the case of photogravitational restricted three body 

1/3 

problem we have ri^ = qi = S and = 1. We 
suppose that due to P-R drag and oblateness, pertur- 
bations in ri(j, are ei, £2 respectively, where e/s are 
very small quantities, then 



gy'(l + ei), 



f2 



1 + 62 



(14) 



Putting these values in equation ( [T^ and neglecting 
higher order terms of small quantities, we get 



- 1 + 362 



3 . 15 ■ 
-A2 + Y^2e2 



362 ( 1 + 



My 



nWi, 



Fig. 10. — The Xi coordinate of equilibrium point when 
0<A2<l,0<gi<l&At = 0.00003, black dot on 
figure indicates the value of X4 in classical case 




Fig. 11. — The j/4 coordinate of equilibrium point when 
0<A2<l,0<gi<l&^ = 0.00003, black dot on 
figure indicates the value of 2/4 in classical case 



£2 



tVKi(l-|A2) 

3/iyo 
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Table 2 x co-ordinate of L4 



A2 


Xi'. qi = l 


: qi = 0.75 


X4 : qi = 0.5 


X4 : qi = 0.25 


X4 : qi = 


0.0 


0.499949 


0.412692 


0.314934 


0.198383 


Complex Infinity 


0.25 


0.374929 


0.30949 


0.236175 


0.148765 


Complex Infinity 


0.5 


0.249906 


0.206284 


0.157412 


0.099146 


Complex Infinity 


0.75 


0.12488 


0.103076 


0.0786479 


0.0495247 


Complex Infinity 


1 


-0.000148407 


-0.000134579 


-0.000118667 


-0.0000980564 


Complex Infinity 



Table 3 


y co-ordinate of L4 








A2 2/4 : gi = 1 


yi-. qi= 0.75 


y4 : gi = 0.5 


^4:91= 0.25 


2/4 : gi = 


0.0 


0.86605 


0.809425 


0.728552 


0.597927 


Indeterminate 


0.25 


0.847795 


0.790465 


0.709777 


0.581035 


Indeterminate 


0.5 


.790543 


0.73068 


0.650249 


0.527174 


Indeterminate 


0.75 


0.684579 


0.61834 


0.536576 


0.422435 


Indeterminate 


1 


0.499831 


0.412542 


0.314727 


0.198036 


Indeterminate 



where j/o is y-coordinate of triangular equilibrium 
points in photogravitational restricted three body prob- 
lem case. Using values ri, r2 with equation ( [T2|) . we 
get following: 



V —y^ + N 



or 



Xi = Xo 



2 Xo 



nH^i[(l-M)(l + f^2)+M(l-#)f] 



3/x(l - ^i)yaxo 



(15) 



- e2)(l + 2 ' 

nW^{rl-rl + l) _^ 

[n" - (1 - /i) + 3(1 - - M 



+ 3e2 - - 562)2^2 



nW^ 



we obtained ei 



nWi 



6(1 - /i)j/o 



— ^ 4z{ r2 = 1 + ^(1- ^A2> 

6(l-M)yo 2 (' ^ 3^11/0^ 2 ^/ 



Now we have x -\- fj, 
therefore, 



(x + fi) 



1 



1 + 



nWi 



3Myo I + M2 



= rf - (a: 

nW^i A2 
6(1 - f^)yo 2 



+ 1 



or 



1- 



6(1 - fi)yo 

nI^i[(l-M)(l + |A2)+Ml-t)f] 
3^(1 - n)yo 



2/4 = 2/0 



S\l-^)A2 



34.)| + 7(l-,.)^lnV2 



3pt(l - ^j)yl 



16) 



where (xo,2/o) are coordinates of L4,L^ in photogravi- 
tational restricted three body problem case as: 



X 1/3 
6 = g/ 



The position of ^4(5) is given by equations (EH), f fTCj) 
which are valid for W"i <C 1 , A2 <C 1 . For simplicity we 
suppose 7=1 — 2/x, gi = 1 — e, with |e| << 1 then the 
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coordinates {x^, 1/4) of ^4(5) can be written as follows: 



X4 



1 
2 



(9 + 7) 



A2 

— + 
2 

nWi 



A2e 
3 



476 



nWi 



6V3 27%/3 

VSf 2e A2 2A2e 

2 I 9 3 9 
(1 + 7)„TT. 47e 



(17) 



Figures ( [IHl E]) and tables ( E]) show how the 
coordinates (x4,y4) points are decreasing functions of 
A2,qi,Wi. Black dots on the figures indicate the 

(0:4,2/4) = (l-Ai,^)• 



4 Comments on the Linear Stability 

In order to study the linear stability of any Lagrangian 
equilibrium point Li[i = 1 — 5) the origin of the co- 
ordinate system to its position [xi,yi) by means of 
X = Xi + a, y = Ui + P, where a = £,e^* , P = rje^* are 
the small displacements f , ri, X these parameters, have 
to be determined. Therefore the equations of perturbed 
motion corresponding to the system of equations ( [3]) , 
(HI may be written as follows: 

a - 2n/3 = aW,, + pUi^y + aW,, + /jf/^ (19) 
i3 + 2na = aUl,+f3Uly + AUl^+pUly (20) 

where superfix i is corresponding to the Li{i = 1 — 5) 

(A^ - XUl^ - Ul,)^ + [-{2n + Uly)X - Uj,y]rj = 



(21) 



[(2n - ui,)x - c/;je + (A^ - xu;^ Wjr^ = 



(22) 



Now above system has singular solution if, 



A2 - XUl 



Ul 



-(2n+t/^.)A-[/i, 



(2n-[/^)A-c/^, x^-xu;^-u; 



At the equilibrium points equations ([3]),(|4]) gives us 
the following: 



b = 2n - - 



3^2 2Wl 



' li 



'2i 



-a(l + e) 



,2 ^5 



rr -r: 



IV 2i 



5^2 



+ 2/, - 



3m^2 



I 9A^(1-M)gi _^ 5^2 



12 \ 2 

2i 



1 2 2t 



^1 , 
^3 "I" ^3 



The points Li,L2,^3 no longer lie along the line 
joining the primaries, since the condition is not sat- 
isfied for them, so taking y 0, ^ — > because 

1/3 

y » Wi,x » Wi, from ( [T]) we have ri « 
In this case fi > 1 for i = 1,2,3, so that for each 
Li,L2, L3 characteristic equation ( [23]) has at least one 
positive root, this implies that these points are unstable 
in the sense of Lyapunov. 

4.1 Linear Stabihty without P-R Effect 

Now we consider the problem when P-R effect is not 
included(VFi = 0), then n = ^^^(1 - ^),r2 = 1. 
The coordinates of triangular points ^4(5) are as: 



' 1/3 



1/3 

^1 



1- 



2/3' 



1/3 A 
gl A2 

2 

1/2 



(24) 



2/3 ' 



1 + 



2/3' 
9i 



Ao 



(25) 



A'' + aX^ + bX^ + cX + d = 



(23) 



a = 0,c= 0,, 



,(* = 4,5) 



b = — 311A2, d = 9/i(l — 



where g = {1 — A2 



1 



istic equation ( [^5]) we obtained, 
A^ 



._-b±{b^-Adf' 



From character- 



(26) 
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For stable motion < 4d < b^, i.e. 
(n^-SMa)^ > 36^(1-^^)5 

In classical case A2 = 0, qi = 1, Wi =0, n = 1, we 
have following: 1 > 27/i(l - m) ^ M < 0.0385201. 
The possible roots of equation ( [26|l are given in ta- 
ble ([H]), we see that all the roots are purely imaginary 
quantities. Hence the triangular equilibrium points 
are stable in the sense of Lyapunov stability provided 
< mouth = 0.0385201. 

Equation ( I26p has imaginary roots Ai^2 = ii'j-'i, 
A3, 4 = ±iuJ2, i = V— 1, this gives us: 



-b ± _ Ad) 



1/2- 



1/2 



(27) 




0.01 0.015 0.02 0.025 0.03 0.035 0.04 

Fig. 12. — The linear stability region for ji — qi pa- 
rameter space and resonance curves uji — kuj2 = 0, fc = 
1,2,3. at A2 = 0,0.02(doted lines) & /i = 0.00003, 
< gi < 1 



There are three main cases of resonances as: 



uji — kuj2 = 0, fc = 1, 2, 3 



(28) 4.2 Linear Stability with P-R Effect 



For fc = 1 we have positive stable resonance and for 
fc = 2,3 we have unstable resonances. Using {^7} and 
( UHl) we obtained a root of mass parameter. 



35 + 2KA2 + SKAj 



^^9g~4K + 9KAl 



Now consider the problem when P-R Effect is included 
i.e., qi ^ l,Wi 7^ 0, A2 7^ 0. Using Ferrari's theorem 
the roots of characteristic equation ( [23|) are given by: 



6{g + KAl) 



(29) 



where K 



(fc2^]^)2 ■ Now we suppose qi = 1 — e, with 
|e| << 1, neglecting higher order terms, we obtained 
the critical mass parameter values corresponding to fc = 
1, 2, 3 as : 

Ati = 0.0385208965 + 0.6755841373^2 

-0.0089174706e (30) 

H2 = 0.0242938971 + 0.4322031625^2 

-0.0055364958e (31) 

H3 = 0.0135160160 + 0.2430452832^2 

-0.0030452832e (32) 

The linear stability region and corresponding main 
resonance curves in fi — qi parameter space are shown 
in figure ( I12p the doted lines are corresponding to 
A2 = 0.02, the curve corresponding to fc = 1, [qi = 
1,^2 = 0,^1 = URouth = 0.038521) is actual 
bounda ry of the stability region these resu l ts are sim- 
ilar to Markellos. Papadakis. and PerdiosI ( 1996[ l and 
others. The critical values of mass parameter fi 
are given in the tables ( HI [5|) for various values of 
qi, ^2- The classical critical val u es of /x are similar 



Deprit and Deprit-Bartholomj (1967). We observe 



to 

that the effect of radiation pressure reduces the linear 
stability zones, these are also affected by oblateness of 
second primary. 



(33) 



where A = ±y/8l-Ab + a^, B = 1{1 + - ^, i 
1,2,3,4 and / is any real root of the equation 



8^^ - Abf + (2c 



This can be written as: 



)l + d{Ab - a"^) - = 



2;2 _ 5/2 _2dl + db=^{{l + ef - 2(1 + e)l + d} 

(34) 



This equation has an exact real root / = ^ for a = 0. 



When a 7^ the roots of characteristic equation ( | 
will be obtained in form of the rapidly convergent series 



2i 



(35) 



Using ( [M 
get ^ = a ± VI + Sai, 



taking the coefficients of only, we 



+ aiaA (l± v/l+8ai) T 



B 



(l+e)(l + e2-6) + d _ 
"^ = ^ 2ib^-4d) 



1 + e 



(36) 
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Table 4 A2 







k ■ qi = l Mfe : 91 = 0.75 ^ik : qi = 0.5 ^fe : 91 = 0.25 /i^ : gi = 



1 0.0385209 

2 0.0242939 

3 0.013516 

4 0.00827037 

5 0.0055092 



0.0363201 
0.0229262 
0.0127632 
0.0078121 
0.00520474 



0.0341355 
0.0215661 
0.0120136 
0.00735548 
0.004901287 



0.0318518 0.0285955 

0.0201415 0.0181056 

0.0112275 0.0101021 

0.00687629 0.00618979 

0.00490128 0.00412616 



Table 5 = 0.02 



k 


/i/c : qi = l 


fik ■ qi = 0.75 


^fe : gi = 0.5 


^ik : qi = 0.25 


M/c : 91 = 


1 


0.0413469 


0.0390477 


0.0367581 


0.0343567 


0.0309186 


2 


0.026094 


0.0246633 


0.0232361 


0.0217366 


0.019585 


3 


0.0145252 


0.0137369 


0.0129496 


0.0121214 


0.0109312 


4 


0.00889015 


0.00841007 


0.00793029 


0.00742525 


0.00669897 


5 


0.00592287 


0.00560384 


0.00528491 


0.00494909 


0.00446598 



We obtained the characteristic roots 



Ai,; 



1 + Vl+8ai 



^3,4 



(37) 

a (1 - Vl + «ai) 



; „Ml-vTT^) .^^ (38, 



where 



1+e 



1+e 



Using above equations ( 157 1 155 ]) . we obtained the roots 
of characteristic equation ( I23p which are presented in 
table ( [7|) for various values of 91 , ^2 ■ We see that at 
least one of the roots Xi{i = 1,2, 3,4) have a pos i- 
tive real part due to P-R effect as in Chernikov ( 1970l ). 
Hence the triangular equilibrium points are unstable in 
the sense of Lyapunov stability. 



motion around the Sun but it is allowed to become a 
satellite of the Earth. This does not mean that it will 
become permanently a satellite of the Earth since its Ja- 
cobi's constant is such that it is not confined to the zero 
velocity oval around Earth. The particle with C > C2 
then it can not change position from the vicinity of Sun 
to that of Earth. The particles with very high C val- 
ues have low relative energy levels and they either move 
around one of the primaries or move for outside of the 
system. Finally we have conclude that, if one the pri- 
mary is exerts the light radiation pressure and second 
primary is an oblate spheroid, then the gravitational 
radiation force and oblateness influence the existence, 
location of equilibrium points. In classical case when 
qi, A2 = 0, we have the coUinear points [Li, L2, L3) 
and two triangular equilibrium points L4, L^. We have 
seen that Li, L2, are linearly unstable, while L4, 
are conditionally stable in the sense of Lyapunov when 
P-R effect is not considered. The effect of radiation 
pressure reduces the linear stability zones while P-R 
effect induces an instability in the sense of Lyapunov. 
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5 Conclusion 



The distances of Li,L2 from the second primary are 
monotonically increasing with fi . The Jacobi's con- 
stant at the iijLs increases monotonically with /i. If 
C < Ci, the particles can leave the system while if 
Ci < C < C2, i.e. the third body is not confined to its 
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Table 7 Roots of characteristic equation when P-R effect is considered 





-^1,2 


: As = 


^3, 


4 : Aa = 


Ai,2 : Aa = .02 


'^3,4 


A2 = .02 


Ai,2 : A2 = .04 


-^3,4 


A2 = .04 


1 




0.0 ± 0.0242i 




0.0 ± 0.9998i 


±0.2522 




0.0 ± 1.0455i 


±0.3692 




0.0 ± 1.0926i 


0.8 


-1.1614 X 


10"^ ± 0.0254i 


2.5627 X 


10"^^ ± 0.9997i 


±0.2523 


-3.4918 X 


10"" ± 1.0455i 


±0.3696 


-7.0611 X 


10"" ± 1.0927i 


0.6 


-2.8140 X 


10"^ ± 0.0272i 


6.9284 X 


10"^^ ± 0.9996i 


±0.2523 


-8.4597 X 


10"" ± 1.0455i 


±0.36700 


-1.7126 X 


10""' ± 1.0928i 


0.4 


-5.5311 X 


10"^ ± 0.0298i 


1.5839 X 


10"^^ ± 0.9996i 


±0.2522 


-1.6621 X 


10"^" ± 1.0454i 


±.3704 


-3.3695 X 


10""' ± 1.09291 


0.2 


-1.1707 X 


10"** ± 0.0347i 


4.3264 X 


10"^^ ± 0.9995i 


±0.2519 


-3.5125 X 


10"^" ± 1.0454i 


±0.3707 


-7.1369 X 


10"^° ± 1.09301 


0.0 




Indeterminate 




Indeterminate 


± 0.0242i 




Indeterminate 


Indeterminate 




Indeterminate 



Note. — 



Tablc[7]is presents the roots of characteristic cguationf I23| l for fi — 0.00003 when P-R effect is included. 



